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Abstract 

In a nanographene ring with zigzag edges, the spin-polarized state and the charge-polarized state are stabilized by the 
on-site and the nearest neighbor Coulomb repulsions, U and V , respectively, within the extended Hubbard model under the 
mean field approximation. In a Mobius strip of the nanographene with a zigzag edge, U stabilizes two magnetic states, the 
domain wall state and the helical state. Both states have ferrimagnetic spins localized along the zigzag edge while the former 
connects the opposite ferrimagnetic orders resulting in a magnetic frustration forced by the topology and the latter rotates 
the ferrimagnetic spins uniformly to circumvent the frustration. The helical state is lower in energy than the domain wall 
state. On the other hand, V stabilizes another domain wall state connecting the opposite charge orders. 
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1. Introduction 

Nanoscale materials are strongly affected by their 
geometries. Carbon nanotubes behave as semicon- 
ductor or metal depending on the diameter and the 
chirality. They opened a new field in fundamental 
science and have become important materials for nan- 
otechnology devices[l,2,3]. Nanographene ribbons, 
nanoscale graphite ribbons, also show strikingly dif- 
ferent characters depending on the shapes of their 
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peripheral edges[4,5]. The nanographene ribbon with 
zigzag edges (zigzag ribbon) has the pecuhar non- 
bonding molecular orbitals localized mainly along the 
zigzag edges (edge states), which cause a sharp peak 
in the density of states at the Fermi energy, while 
such states are completely absent in the ribbon with 
armchair edges. Moreover, zigzag ribbons will show 
rich natures due to the competing spin and charge 
orders: the introduction of the on-site Coulomb repul- 
sion U under the mean field approximation induces 
the Fermi-instability of the edge states and results in 
the spin-polarized (SP) states with the opposite fer- 
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rimagnetic orders localized along the zigzag edges [4]; 

On the other hand, the nearest neighbor Coulomb 
repulsion V induces another Fermi-instability and re- 
sults in a novel ferroelectrically charge-polarized (CP) 
state with the unlike charges localized along the zigzag 
edges and competes with the SP state [6]. Recently, 
microscopic Mobius strips of NbSes were synthesized 
[7] and stimulated theoretical studics[5,8,9]. 

A Mobius strip has unique topological natures dif- 
ferent from a ring; the former has only one side and 
one edge while the latter has two sides and two edges. 
Then, how such an exotic topology affects the elec- 
tronic states in nanographene? This is the question 
that wc study here. 

In this paper wc demonstrate that the Mobius strip 
has both the SP domain wall (SP-DW) state [5] and the 
helical state due to U. We also show that V stabilizes 
the novel domain wall state connecting the opposite 
charge orders of the CP state (CP-DW). 



2. Model and Methods 

Figure 1 illustrates the geometry of a zigzag ribbon 
with an inversion symmetry to make a Mobius strip or 
a ring. N and L are the width and length of the ribbon. 
Wc consider the case both A'^ and L to be even. Note 
that the twofold coordinated sites locate only at the 
zigzag edges, PR and QS, where the edge states are 
mainly localized. Since the zigzag ribbon is bipartite, 
A and B sites are assigned by black and white circles, 
respectively. The Mobius strip is obtained by twisting 
one end of the ribbon, e.g., PQ and joining the ends so 
that the P site is connected with the S site and the Q 
site is connected with the R site, defining the Mobius 
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Fig. 1 . The schematic structure of a bipartite nanographene 
ribbon with zigzag edges to make a non-bipartite Mobius 
strip or a bipartite ring; • (o), A (B) site in the bipartite 
lattice. Dashed rectangle denotes a unit cell in a ring. 



boundary condition. 

Note that both the P and S sites belong to A sub- 
lattice while both the Q and R sites belong to B sub- 
lattice. Since the Mobius boundary condition forces to 
connect the same sublattice, the Mobius strip is not 
bipartite, while the periodic boundary condition keeps 
the system bipartite. 

We treat a half-filled tt electron system on the zigzag 
ribbon with the extended Hubbard Hamiltoriiari, 

Heh = - * X] ('^L'^J'.'' + 

1 1 

+ C/^(ni,T--)(ni,i--) 

i 

+ l/^(ni-l)(n,-l). (1) 

(id) 

Here, (i, j) denotes a pair of neighboring carbon sites 
and —t is the transfer integral between them, cj ^ cre- 
ates a TT electron with spin s on the i-th site and Ci,s 
annihilates the one. is the corresponding number 
operator, and rii = ni,s- 

To study the possibility of the helical state, we gen- 
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Fig. 2. Spin density distributions of (a) the spin-polarizcd 
(SP) state in a ring, (b) the SP domain wall (SP-DW) state 
and (c) helical states in a Mobius strip with iV X L = 4 X 20 
sites for (7 = f and V = 0. Arrows represent spins of tt 
electrons. In (a) and (b), only z component of spin densities 
Si^z is finite: In (c), y component st^y is zero. 
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Fig. 3. U dependences of the MO energies of the SP (solid 
hne), SP-DW (Q) and heUcal (•) states for V = 0. 

eralized the unrestricted Hartree-Fock approximation 
to set a two body operator as [10], 

- (cl^Ci,-j){cl^Ci,i) +{cl^Ci,i){cl^Ci,^). (2) 

The u{= X, y, z) component of the spin density and 
the charge density at i-th site are given by Si,„ = 
^a.i3 (4,a('^i')a/3Ci,^i>/2anddi = 1 -(ni), respectively, 
where CTu is the Pauli matrix. 



3. Results and Discussions 

Figure 2 shows spin density distributions of (a) the 
SP state in a ring, and the (b) SP-DW and (c) heli- 
cal states in a Mobius strip with NxL = 4x20 car- 
bon sites for [7 = t and F = 0. In the SP state, the 
up- and down-spin ferrimagnctic orders appear along 
the upper and lower zigzag edges, respectively, because 
the spin-density-wave like spin-orientation order is sta- 
bilized by U in the bipartite system with the edge 
states at twofold sites [4]. The Mobius boundary con- 
dition forces to connect the opposite ferrimagnctic or- 
ders along the zigzag edge and causes the magnetic 
frustration. It allows the SP-DW state with the op- 
posite ferrimagnctic orders connected by the interface 
where the magnitudes of spins are suppressed to re- 
duce the frustration [5]. On the other hand, the helical 
state keeps the magnitudes of spins along the zigzag 
edge with the ferrimagnetic feature in a unit cell. In 
this state, the directions of spins at twofold sites ro- 
tate uniformly along the zigzag edge to circumvent the 
magnetic frustration due to the Mobius boundary con- 
dition. The spin rotates by 2tt in the SxSz plane dur- 



h'ig. 1. (a) U dependence of the energy difference 
AE = -Esp-DW ~ ^helical ^OT the System with 4x40 sites 
and V = 0, where Esp-dw and ^helical are the energies of 
the SP-DW and helical states, respectively, (b) Dependence 
of AE per unit cell on the ribbon length L for U = t,V = 
and N = A. 



ing the travel along the zigzag edge to return to the 
starting point, showing the helical ferrimagnetic or- 
der. Therefore, the SP-DW and helical states satisfy 
the constraint of the Mobius boundary condition con- 
necting the opposite zigzag edges of the ribbon, i.e., 
the opposite ferrimagnetic domains in the SP state. A 
transformation of all spins in the helical state under 
an Eulcr angle rotation in the spin space causes an- 
other degenerate helical state, since the spin space is 
isotropic in this model. Total magnetization are zero 
in the SP-DW and helical states. 

In the case of much larger U (~10t), the magnitudes 
of spin densities become almost the same at all sites in 
the SP and helical states as well as in each ferrimagnetic 
domain of the SP-DW state: the ferrimagnetic feature 
turns into the antiferromagnetic one with increasing U . 

Figure 3 shows the U dependences of the molecular 
orbital (MO) energies around the Fermi level in the SP, 
SP-DW and helical states for the system with NxL = 
4x40 sites. In the SP state, U opens the energy gap be- 
tween the degenerate edge states at the Fermi energy, 
while in the SP-DW state the frontier orbitals remain 
inside the gap [5] similar to the soliton state in poly- 
acetylene [11]. On the other hand, in the helical state 
the frontier orbitals behave more similarly to that in 
the SP state than in the SP-DW state. This is because 
the variation of the spin direction in the helical state is 
extended over the system. In the shorter ribbon with 
NxL = 4x20, the difference between their behaviors 
in the SP-DW and helical states is not so appreciable 
because a half pitch L in the helical state is closer to 
the size of the interface in the SP-DW state. 

For U > 0, the helical state is lower in energy than 
the SP-DW state, as shown in Fig. 4 (a). This is consis- 
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Fig. 5. Charge density, di, distribution of (a) the 
charge-polarized (CP) state in a zigzag ring and (b) the do- 
main wall state connecting the opposite charge orders of the 
CP state in a Mobiiis strip for U = 0.2t and V = OAt; •(o), 
positive (ncgati\'c) cliarge densities. The radius of each cir- 
cle denotes tlic magnitude of the charge density; the max- 
imal ones arc 0.23 in (a) and 0.24 in (b). The system sizes 
are NxL = 4x20 sites. 

tent with the MO's behavior shown in Fig. 3: the wider 
gap and the enhanced decrease in the occupied MO en- 
ergies in the helical states than in the SP-DW states. 
Figure 4 (b) shows that the energy difference per unit 
cell is smaller in the shorter ribbon, as L approaches 
the size of the interface in the SP-DW state. 

We recently reported that V stabilizes the novel CP 
state with a finite electric dipole moment under pe- 
riodic boundary condition[6]. Figure 5 (a) shows the 
charge density distribution of this state for U = 0.2t 
and V = OAt. The positive and negative charges are 
localized along the lower and upper zigzag edges, re- 
spectively, due to V and they form a finite electric 
dipole moment pointing from the upper edge to the 
lower edge. Net charge along each zigzag edge arises 
from much larger charge density at twofold sites than 
at threefold sites due to the edge states. The magnitude 
of the electric dipole moment per unit cell is 6.2 Debye 
in the parameters. This exotic electric dipole moment 
is induced by V and breaks the inversion symmetry of 
the system. The CP state competes with the SP state. 

In a Mobius strip, corresponding domain wall state, 
CP-DW state, is stabilized as shown in Fig. 5 (b). It 
connects the opposite charge orders of the CP state. 
The Mobius boundary condition allows only the CP- 
DW state in this case since charge is a scalar, while 
it allows the SP-DW and helical states in the spin or- 
dered case since spin is a vector. The phase diagram 
for their competition in the U-V space will be reported 



elsewhere. 



4. Conclusion 

In summary, we demonstrated that a Mobius strip of 
nanographene with a zigzag edge has the spin-polarized 
domain wall state and the helical ferrimagnetic state 
caused by the on-site Coulomb repulsion within the 
extended-Hubbard model under the mean field approx- 
imation and the former is lower in energy than the lat- 
ter. It also has the charge-polarized domain wall state 
stabilized by the nearest neighbor Coulomb repulsion 
which connects the opposite charge orders. 
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